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ABSTRACT 


The  strain  energy  of  a  spheroidal  inclusion  was  evaluated 
exactly  using  the  Eshelby  theory.  Numerical  results  for  an 
oblate  spheroid  are  presented  in  a  parametric  form  in  terns  of 
the  transformation  strain  tensor.  Using  atomisitc  transformation 
mechanisms,  the  transformation  strain  was  determined  for  bcc-hcp 
martensitic  transformation  in  Ti  and  its  alloys.  The  lattice 
correspondence  satisfied  the  Burgers  relationship  and  the  c/a-natio 
in  the  preduct  phase  was  taken  as  1.586.  The  habit  plane  was  pre¬ 
dicted  on  the  basis  of  the  strain  energy  minimization  principle. 
Results  of  the  calculation  indicate  that  the  strain  energy  is 
minimized  when  the  morphology  of  hep  Ti  martensite  is  a  thin 
disc-shaped  inclusion  lying  on  a  plane  close  to  (9  9  13^,  in  ^ 

excellent  agreement  with  experimental  observations. 

The  present  approach  is  discussed  in  detail  and  compared  with 
the  crystallographic  theory  based  on  the  invariant  plane  strain 
hypothesis. 


i.  iMnmmoN 


Because  of  the  technological  importance  of  ferrous  martensites, 

various  aspects  of  nartensitic  phase  transformation  have  been  investi- 
fl  2) 

gated  extensively.  *  The  phenomenological  crystallographic  theory 
of  martensite  formation  was  developed  in  Idle  1950's  and  is  based  on  the 
concept  of  the  invariant  habit  plane.  The  theory  resolves  the  total 
shape  deformation,  F,  into  a  pure  strain  that  converts  the  parent  into 
the  product  lattice,  P,  a  lattice  invariant  shear  "that  produces  two 
undistorted  planes,  5,  and  a  rotation  that  ensures  the  invariant  inter¬ 
face  plane,  R.  The  basic  hypothesis  of  the  theory  is  that  the  total 
shape  deformation  is  an  invariant  plane  strain.  The  whole  of  the  volume 
change  of  the  transformation  is  produced  as  an  expansion  or  contraction 

nomal  to  the  habit  plane  as  well  as  a  small  dilatation  within  the  habit 

(3) 

plane  in  the  Bowles-ffackenzie  (B-M)  formulation.  The  various  formila- 

(3) 

tions  of  the  theory  were  presented  by  Bowles  and  Mackenzie,  Bullough 

and  Bilby,^and  Wechsler  et  al.^  and  in  a  modified  form  by  Ross  and 
(6) 

Crocker,  et  al.  These  formulations  lave  been  successful  in 
accounting  for  the  crystallographic  characteristics  of  nary  transformat  ions , 
such  as  ferrous  transformations  with  the  {3  15  10}p  habit  plane  and  the 
bcc-hcp  transformation  in  Ti  and  Zr  with  the  {334 }fi  habit  plane. ^ 

However,  a  variety  of  martensites  are  difficult  to  reconcile  with  the 
theory.  This  group  includes  ferrous  martensites  with  the  habit  plane  of 
{225}p,  {259}p,  {142 }p,  etc.  and  Ti  martensites  with  the  {344 }g  habit 
plane. it  is  also  rat  clear  from  the  theory  vhich  factors  are 
responsible  for  various  crystallographic  characteristics  observed  in 
different  alloys. 


2 


It  is  well  known  that  the  elastic  defamation  in  the  particle  of 

new  phase  and  in  the  parent  lattice  controls,  to  a  large  extent,  a 

martensitic  transformation  and  the  coherent  stage  of  a  phase  trans- 
(8) 

formation.  In  the  above-mentioned  crystallographic  theory,  it  is 
assumed  that  a  thin  plate  martensite  with  the  undistorted  interfaces 
satisfies  the  criterion  of  a  strain  energy  minimum.  This  assumption 
is  valid  in  the  cases  of  plate  martensites.  Christian  considered 
accommodation  strains  in  the  fomation  of  a  martensite  plate  which  is  an 
oblate  spheroid.  The  treatment  is  based  on  the  Eshelby  theory  of  the 
strain  energy  in  a  constrained  transformation/10^  but  is  valid  only 
for  a  large  diameter-to-thickness  ratio.  His  result  shows  that  the 
strain  energy  per  unit  volume  tends  to  zero  with  an  increasing 
diameter-to-thickness  ratio.  It  is  clear  that  more  exact  evaluations 
of  the  strain  energy  via  a  generalized  elasticity  treatment  are 
necessary  in  order  to  develop  a  new  theory  that  is  based  on  the  principle 
of  minimum  free  energy.  FUrthennore,  the  atomistic  mechanisms  involving 
dislocations,  stacking  faults  and  twins  must  be  integrated  into  the  new 
theory  so  that  only  the  probable  modes  of  deformation  are  to  be  considered. 

In  this  paper,  we  establish  the  correspondence  between  the  crystal¬ 
lographic  and  Eshelby  theories  and  describe  the  calculation  of  the  strain 
energy  of  spheroidal  inclusions  using  the  Eshelby  theory.  The  results 
are  combined  with  the  consideration  of  atomistic  mechanisms  in  a  detailed 
analysis  of  martensitic  transformation  in  Ti  alloys.  It  is  shown  that 
the  results  of  strain  energy  minimum  determination  correlate  well  with 
experimental  observations  on  the  habit  plane.  Other  features  of  the  present 
approach  is  discussed. 


2. 


Theoretical  Considerations 
2.1  The  Eshelbv  Theory 

Eshelby^10^  has  obtained  a  general  solution  for  the  elastic  state 
of  an  ellipsoidal  inclusion  and  the  matrix,  when  the  inclusion  within 
an  isotropic  elastic  solid  undergoes  a  change  of  shape  and  size,  which 
would  be  homogeneous,  but  for  the  contraint  imposed  by  the  matrix.  Be¬ 
cause  of  the  presence  of  the  matrix,  internal  stresses  exist  both  inside 
and  outside  the  inclusion.  In  order  to  find  the  elastic  field,  Eshelby 
uses  a  set  of  five  imaginary  operations,  as  shewn  in  Fig.  1.  These 
are 

(a)  Cut  around  the  inclusion  which  is  to  transform  and  remove 
it  from  the  matrix, 

(b)  Allow  the  transformation  to  proceed  without  the  constraint; 
namely,  stress-free  transformation  defined  by  e^j, 

(c)  Apply  surface  tractions  to  restore  the  original  shape  of  the 
inclusion, 

(d)  Replace  the  inclusion  into  the  hole  in  the  matrix  and  rejoin  the 
material  across  the  cut, 

Ce)  Spring  back  to  the  self  stress  state  by  the  release  of  the 
surface  tractions. 

The  net  effect  of  these  operations  is  essentially  equivalent  to  the  total 
shape  deformation,  F,  of  the  crystallographic  theory;  that  is,  Operation  b 
corresponds  to  the  lattice  deformation,  RP,  and  Operation  c  to  the  lattice 
invariant  shear,  3.  Since  the  crystallographic  theory  deals  with  the 
constrained  condition,  Operations  a  and  d  are  omitted.  Rjrther  the 
assumption  of  a  very  thin  plate  shape  leads  to  vanishing  strain  energy 
so  that  the  relaxation  process  of  Operation  e  needs  not  be  considered. 


On  the  other  hand,  only  pure  shear  strain  (or  its  combination)  is 
allowed  for  §  in  contrast  to  any  homogeneous  strain  permissible  in 
Operations  b  and  c.  This  leads  to  the  main  distinction  between  the  two 
approaches.  The  Eshelby  theory  provides  no  shape  change  after  Operations 
d  and  a  limited  shape  change  or  rotation  following  Operation  e,  whereas 
the  total  shape  deformation ,  F  is  merely  reduced  to  the  state  of  an  in¬ 
variant  plane  strain  in  the  crystallographic  theory. 

The  total  strain  energy  in  matrix  and  inclusion,  E,  is  obtained 


total  vol 


°ij  eij  dV=  "  ij  eij  ^ 

•  v  T 


where  oA  is  the  stress  in  the  inclusion,  and  Vj  is  the  volume  of  the 
inclusion,  and  the  summation  of  repeated  indices  is  implied.  In  the 
case  of  an  ellipsoidal,  inclusion,  cn?  in  it  is  uniform  and  is  given  by 


rij  =  Cijk&(Sk2jnn  emn  "  ekl) 


where  the  Ch..^  are  the  elastic  stiffness  and  the  are  the  coefficients 

relating  the  transformation  strain  and  the  constrained  strain  by 

ekJl  =  SkAmn  emn  *  (3) 

While  the  evaluation  of  E  is  cumbersome,  the  solution  to  the  elasticity 
problem  of  transformation  exists  in  principle.  Thus,  the  remaining  tasks 
are  to  obtain  suitable  e.t  for  a  specific  transformation  mechanism  and  to 
establish  a  strain  energy  and  free  energy  mininum  condition(s)  which  depends 
strongly  on  geometry. 


in  terms  of 


2.2  The  Strain  Energy  of  Spheroidal  Inclusions 

Eshelby  described  the  general  solution  for 
elliptic  intergrals,  from  which  E  can  be  obtained.  To  date,  however, 

E's  for  only  few  special  cases  have  been  evaluated.  While  the  Eshelby 
theory  is  valid  for  general  ellipsoids,  we  shall  treat  only  spheroidal 
inclusions.  This  greatly  simplifies  the  calculation  without  sacrificing 
the  applicability  to  phase  transformation  problems.  It  is  also  assumed 
that  the  elastic  constants  of  the  inclusion  are  identical  to  those  of  the 
matrix  and  that  only  one  inclusion  is  present  in  an  infinite  media. 

This  section  presents  the  solutions  of  E's  for  several  simple  cases 
and  a  general  solution  for  an  oblate  spheroidal  inclusion.  Consider  a 
spheroidal  inclusion,  the  major  axes  of  vhich  coincide  with  the  coordinate 
axes;  namely 

(x-j/a)2  +  (x2/b)2  +  (x3/c)2  =  1  (4) 

Two  of  the  constants  are  identical,  e.g. ,  a  =  c.  The  aspect  ratio  k  is 
defined  as  k  =  b/a  =  b/c  for  this  case.  Fig.  2  shows  such  an  inclusion 
with  k  <  1  (an  oblate  spheroid) .  In  the  expressions  below,  Poisson's  ratio 
is  always  taken  to  be  1/3. 

Case  1. 

"t  t  It 

When  the  only  nonzero  components  of  are  e^2  =  ^i*  s-traih 
energy  E^  is  given  by 

E^  =  cyiCe^Vj,  Vj  =  j  ir  abc  (5) 

where  is  a  constant  that  depends  on  the  orientation  and  the  aspect 

ratio  of  the  inclusion  and  p  is  the  shear  modulus,  can  be  expressed 

exactly  in  terms  of  elliptic  integrals,  but  approximate  expressions  for 


are  listed  in  Table  I.  The  range  of  k  where  these  are  valid  to  with¬ 
in  5%  of  the  exact  values  is  also  given  on  the  basis  of  a  comparison 
with  exact  calculations. 

Case  2. 

When  all  the  components  of  eth  except  e22  vanish,  the  strain  energy 


E2  is  given  by 


E2  =  a2  v(e22)  Vj 


(6) 


The  approximate  expressions  for  a2  are  tabulated  in  Table  I.  E-^  and  E2 
have  no  interaction  term  so  that  the  sum  of  the  two  gives  the  total  strain 
energy  correctly. 

Case  3. 

When  all  the  elements  of  e.^  are  nonzero,  the  strain  energy  E_  is 

1J  O 

given  by  cn?  e^j  Vj/2.  When  an  oblate,  spheroidal  inclusion  with  a  =  c  >  b 
exists  in  the  matrix,  the  strain  energy  is  given  by 


2 

E3  =  U  VI  h(ell  +  e33)  +  a2  ell  e33  +  a3  (ell  +  e33)  e22 
+  a4(e^)  +  a5[(e^)2  +  (e*  J]  +  a6(ej)2)  , 


(7) 


where  the  numerical  values  of  coefficients  a^  through  ag  are  plotted 
against  aspect  ratio  k(=  b/a)  in  Fig.  3.  For  oblate,  spheroids  of 
different  configurations,  e.g. ,  a  =  b  >  c  or  a  <  b  =  c,  Eg  is  obtained 
by  interchange  of  the  indices  of  e^.. 
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2.3  The  Present  Theory — An  Outline 


The  main  objective  of  this  study  is  to  combine  the  atomistic 
mechanism  of  martensitic  transfomation  with  the  generalized  Eshelby 
treatment.  The  present  theory  consists  of  the  following  steps: 


(i)  Determine  the  transformation  strain  (e^)  based  on  the  lattice 
correspondence  and  lattice  parameters  of  the  parent  and  product 
phase  as  well  as  on  the  probable  mode  of  the  et . . 

(ii)  Evaluate  the  total  strain  energy  of  a  spheroidal  inclusion  and 
the  matrix  as  a  function  of  the  shape  and  orientation  of  the 
inclusion  using  the  Eshelby  tensor  obtained  from  eT.  . 

(iii)  Determine  the  strain  energy  minimum  conditions  and  the  habit 

orientations.  The  minimum  free  energy  condition  can  be  deduced 
by  including  the  interfacial  energy  term. 


.  .  t  tr  . 

Although  no  restriction  is  placed  on  e^j ,  generally  e^  is  decomposed 

into  one  or  two  atomic  shear  mechanisns  and  dilatation  terms.  It  is 

assumed  that  both  the  inclusion  and  matrix  are  isotropic  and  homogeneous. 

In  this  paper,  we  describe  a  specific  application  of  the  present  theory 

to  nrartensitic  transfomation  in  Ti  alloys.  Details  of  the  calculation  are 

presented  in  the  following  sections. 


2.4  The  Strain  Energy  Calculation 

When  the  najor  axes  of  a  spheroidal  inclusion  are  parallel  to  the 
coordinate  axes,  E  can  be  evaluated  by  the  procedures  in  Section  2.2, 
once  etj  and  the  inclusion  geometry  are  known.  However,  these  orientations 
are  not  necessarily  the  minimum  energy  configuration.  In  order  to  determine 
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E  for  a  spheroidal  inclusion  of  arbitrary  orientation,  the  coordinate 
axes  are  transformed  to  coincide  the  major  axes  of  the  inclusion.  This 
requires  the  simultaneous  transformation  of  e^  into  the  new  coordinate 
system,  but  enables  the  subsequent  use  of  the  procedures  in  Section  2.2. 

When  e*j  in  the  original  coordinate  system  is  given,  the  rotation 
of  the  coordinate  by  angle  6  about  the  x3  axis,  as  shown  in  Fig.  4, 
results  in  a  new  strain  tensor, 

/cos 0  -sin0  o\  /e ii  ^2  ei3  \  /  cos®  s^ri®  0 

|  sine  cos6  0  jf  e12  e22  e23  j  j  -sine  cos6  0 

\  0  0  1/ \e13  e23  e33/  V  0  0  1  / 


<4j>e  = 


(8) 


Since  the  inclusion  is  spheroidal,  another  rotation  by  angle  ip  about  one 
of  the  new  coordinate  axes,  e.g. ,  the  x^  axis,  produces  the  desired  re¬ 
orientation  effect  and  enables  the  evaluation  of  E.  However,  most  of  the 
components  of  transformation  strain  in  the  doubly  rotated  coordinate 
system  are  nonzero  and,  therefore,  the  most  general  formula  for  E  must 
be  used.  After  the  two  rotation  operations,  a  new  transformation  strain 
tensor  is  given  by 


<eij)6,t|» 


0 

cos ip 
simp 


0 


cosip 

-sirop 


(9) 


3.  Bcc  to  hep  Martensitic  Transformation 

The  transformation  of  high  temperature  bcc  phase  of  Ti  and  Zr  alloys 
during  quenching  is  martensitic  and  results  in  at  least  four  structures, 
including  hep,  fee,  orthorhombic  and  fee  orthorhombic The  hep 
martensite  is  most  common  and  occurs  in  Ti,  Zr,  Ti-Mo  alloys,  Zr-Nb  alloys 
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and  others.  In  pure  metals  and  other  dilute  alloys,  the  lath  martensite 
forms,  whereas  the  plate  martensite  occurs  with  increasing  solute  content . 
The  habit  plane  is  typically  (334)g  or  (8  9  12)g.  The  orientation  relation¬ 
ship  of  the  plate  nartensite  to  the  iratrix  has  been  established  to  be  that 
(12) 

due  to  Burgers,  i.e.. 


( Oil )B//( 0001  >H,  cm]B//[ii2o]H  . 


The  experimental  orientation  relations  in  Ti  and  Zr  also  are  close  to  to 
the  Burgers  relationship. 

The  atomic  movements  for  this  transformation  can  be  accomplished  by 

shearing  on  (211 )g  along  [lll]g  and  a  dilation  along  [2ll]g  as  shown  in 

Fig.  5.  The  stacking  of  atoms  on  (011)B  is  identical  on  every  other  layer 

of  (011)g  and  the  illustrated  atom  movements  result  in  the  correct  atomic 

arrangement  and  stacking  sequence  of  (DOOl)^.  However,  the  atoms  on  the 

middle  layer  must  be  moved  as  well.  This  is  the  so-called  atomic  shuffle, 

and  is  ignored  in  the  present  theory.  This  atomistic  mechanism  was  first 

(12) 

proposed  by  Burgers  and  satisfies  the  Burgers  orientation  relationship. 

When  the  nearest  neighbor  distance  remains  unchanged,  an  ideal  hep  lattice 
is  produced  by  the  atomic  movements  and  the  c/a-natio  is  1.633.  In  the 
subsequent  discussions,  we  shall  always  refer  to  the  directions  [lll]g, 
[2il]g  and  [011]g  as  x^,  x2»  and  xg  axes,  respectively  (cf.  Fig.  2).  These 
correspond  to  the  direction  of  shear,  the  normal  to  the  plane  of  shear  and 
the  normal  to  the  plane  containing  all  the  atomic  movements,  respectively. 
In  terms  of  this  coordinate  system,  the  transformation  strain  for  idle  above 
atomic  movements  is  given  by  -  \ 


.088 

.081 


(10) 


0 


0 


0 
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Since  the  c/a-ratio  of  Ti  nartensite  may  be  approximated  by  that  of  pure 
Ti,  1.586,  additional  dilataticnal  strain  components  must  be  considered. 
Assuming  that  the  interplanar  spacing  of  (011)g  remains  unchanged,  the 
transformation  strain  for  this  case  becomes 


0.030  0.088  0 

0.088  -0.051  0 
0  0  0 


(11) 


The  ratio  of  lattice  parameters  of  the  product  and  parent  lattices, 
ajj/ag,  is  equal  to  0.892  in  this  case,  vfiich  is  slightly  smaller  than  an 
experimentally  observed  value  of  0.899.  The  transformation  strain  of 

(13) 

Eq.  (11)  compares  to  the  one  for  zirconium  obtained  by  Kelly  and  Qroves, 
which  is  given  in  the  present  coordinate  system  as 


'U 


0.033  0.094  0 

,0.094  -0.033  0 

k0  0  0.02y 


(12) 


In  addition  to  e^ ,  the  uniform  rotation  accompanies  the  trans¬ 
formation.  In  the  present  case,  it  is  given  by 


(13) 


The  sum  of  e^  and  is  the  net  distortion  6^j  during  the  trymsfccsation. 
Note  that  an  atomistic  mechanism  of  phase  transformation  defines  8*j , 
as  only  the  e^  term  of  8^  is  considered  in  the  Eshelby  theory. 


where- 


r.mutt^-  '  ‘'>s  /  vm 


Results 


4.1  Simple  Cases 

When  e*j  is  given  by  Eq.  (10)  and  the  coordinate  system  of  Fig.  2 
is  used,  the  strain  energy  is  divided  into  the  shear  component  arid 
the  dilatational  component  Eg.  Using  the  exact  expressions,  E^  and  E2 
have  been  evaluated  and  are  shown  as  a  function  of  k  in  Figs.  6  and 
7,  respectively.  The  orientation  of  an  oblate  spheroid  is  identified 
by  the  normal  to  the  broad  face,  e.g. ,  (2ll)-disc,  whereas  that  of  a 
prolate  spheroid  by  the  long  axis,  e.g.,  [110]-needle.  Figure  6  shows 
that  the  values  of  E^^  are  identical  for  the  (211)-  and  (Ill) -discs  and 
for  the  [211]-  and  [III] -needles.  E^  is  the  lowest  for  these  orientations 
for  disc-like  inclusions  having  snail  k  values,  and  increases  with  increasing 
k  up  to  k  a  2,  after  which  E  decreases  slightly.  The  (Oll)-disc  has  the 
highest  E^,  but  E1  for  the  [Oil]  needles  is  the  lowest  among  -the  needles. 
Figure  7  shows  the  similar  result  for  Eg.  The  lowest  value  is  obtained 
for  the  (3Il)  disc  at  vanishing  k.  Anong  needle-like  inclusions,  the 
[III]  and  [Oil]  orientations  tad  the  lowest  values  of  Eg.  Eg  either 
decreased  or  increased  with  raising  k  for  a  given  orientation.  This  is 
in  contrast  to  the  strain  energy  due  to  dilatation  as  calculated  by 
Nabarro^14^  and  shown  here  by  a  dashed  line.  The  latter  exhibits  a 
maximum  for  the  spherical  shape  and  a  minimum  for  the  thin  discs. 

The  total  strain  energy,  E  =  E^  +  Eg,  is  plotted  against  k  in 
Fig-.  8.  The  lowest  energy  configuration  corresponds  to  a  thin  oblate 
spheroid  with  [!!Il]  normal  to  the  broad  face,  or  the  [5Il]-disc.  The 
second  lowest  energy  configuration  is  the  [Oil] -needle.  The  total  strain 
energy  for  the  [5Il]-disc  continues  to  decrease  with  decreasing  k,  but  that 
for  the  [Oil] -needle  approaches  the  asymptotic  value  at  k  >5. 

While  a  suitable  selection  of  parameters  is  difficult,  the  effect 
of  the  interface  energy  contribution,  E  can  be  assessed  to  a  limited 

S 

*2  2 

degree.  Assuming  that  the  volume  of  an  inclusion  is  S  x  10  cm5 ,  the 
interface  energy  Y  *  20  erg/ an2 ,  and  (i  -  2  x  101  ‘dyn/cnt2 ,  E_  was  evaluated 

5 
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as  a  function  of  k.  The  total  free  energy,  Erp  =  E  +  Eg,  is  shown  in 
Fig.  9.  For  this  snail  inclusion  size,  the  minimum  of  Ej,  exists  at 
k  =  0.37  for  the  (2II)-disc. 

4.2  The  Minimum  Strain  Energy  Configurations 

When  an  oblate  spheroid  with  the  broad  face  nornel  of  [5Il]  (the 
x2  axis)  is  rotated  about  [Oil]  (the  x3  axis)  as  shown  in  Fig.  4;  the 
strain  energy  of  the  spheroid  varies  with  the  angle  or  rotation  0  and 
has  a  two-fold  rotational  symmetry  as  shown  in  Figs.  10  and  11.  In  both 
cases,  several  curves  corresponding  to  different  values  of  k  are  shown. 
Figure  10  was  obtianed  using  the  transformation  strain  of  Eq.  (10),  and 
Fig.  11  using  that  of  Eq.  (11),  respectively. 

Figure  10  represents  the  strain  energy  of  a  spheroidal  inclusion 
transforming  to  an  hep  lattice  with  the  ideal  c/a-natio.  The  difference 
between  the  maximum  and  minimum  values  in  E  tends  to  vanish  with  k  approach¬ 
ing  unity.  The  minor  peak  at  0  =  150°  is  much  smaller  than  the  peak  at 
0  =  55°.  The  minimum  value  of  strain  energy,  and  the  corresponding 

value  of  0  are  summarized  in  Table  II.  Also  shown  in  the  table  are  the 

ranges  of  0  C*A0)  within  vhich  E  differs  less  than  1%  of  the  strain  energy 

t 

for  a  spherical  inclusion  with  the  identical  e^.  The  strain  energy 
vanishes  at  8  *  0°  or  180°  and  115°  for  k  =  0.  When  k  is  raised  to  0.1, 
the  minimum  value  of  the  strain  energy  increases  and  two  strain  energy 
minima  exist  at  0  =  118  and  176.5°.  At  k  *  0.4  and  0.6,  only  one  broad 
minimum  is  found  centering  at  9  *  147°,  and  E  is  nearly  the  same  within 
*  15°.  Obviously,  no  angular  dependence  exists  for  the  spherical  inclusion, 
k  =  1. 
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tabu:  ii 

Strain  Energy  Minina 


1] 

k 

W"vi 

6  (degrees) 

±  A8 

i  Atj>  (degrees) 

0 

0 

115,  180  (or  0) 

2 

4 

Eq. 

CIO) 

0.1 

3.16  x  10‘3 

118,  176.5 

3 

5 

0.4 

8.05  x  10~3 

147 

16 

8 

0.8 

11.9  x  10“3 

147 

15 

15 

0 

0 

9,  106 

2 

4 

Eq. 

(11) 

0.1 

3.05  x  10”3 

8,  106.5 

2 

4 

0.4 

7.64  x  10"3 

4,  111 

6 

7 

0.8 

9.68  x  10“3 

-20,  120 

50 

25 

15 


When  the  transformation  strain  [Eq.  (11)]  corresponding  to  the 


c/a-ratio  of  1.586  in  the  product  lattice  is  used,  the  angular  depen¬ 
dence  of  the  strain  energy  shown  in  Fig.  11  was  obtained.  For  the  k 
values  less  than  unity,  two  strain  energy  minima  always  exist,  and  the 
minor  peak  at  $  =  150°  is  more  prominent  than  in  the  previous  case. 

As  listed  in  Thble  II,  the  minimum  strain  energy  orientations  shifted 
by  9  to  37°  between  the  corresponding  values  determined  using  the  two 


The  values  of  1^^  on  a  single  curve  in  Figs.  10  and  11  are 
identical  to  each  other.  The  variations  of  with  k  are  shown  in 
Fig.  12.  Below  k  —  0.1,  the  values  of  I^n  differ  only  slightly*  but 
the  two  curves  diverge  with  increasing  k.  It  is  also  clear  that  Emi.n 
far  a  given  k  is  lower  when  of  Eq.  (11)  is  employed  to  produce  the 
correct  c/a-ratio  for  titanium  than  when  e^j  of  Eq.  (10)  is  used  for  the 
calculation. 

Results  of  the  strain  energy  claculaticns  for  a  doubly  rotated  in¬ 
clusion  are  sumnarized  in  Figs.  13-15.  For  presentation  of  voluminous 
data,  the  values  of  strain  energy  are  normalized  by  the  maximum  strain 
energy  (E^^)  for  a  given  etj  and  k.  The  magnitude  of  E  is  indicated 
by  intergers  1,  2,  3,  ....  9,  which  correspond  to  (0.10  t  0.01)  ^|L)tv> 

(0.2  *  0.01)  E^,  (0.3  ±  0.01)  1^, - (0.0  ±  0.01)  1^,  rwpec- 

tively.  The  locations  corresponding  to  (l.O  -  O.Ol)  E^  are  indicated 
by  asterik  (*)  marks,  while  that  of  is  given  by  a  plus  (♦)  mark. 
Figs.  13  and  14  show  the  positions  of  various  levels  of  E  as  a  function 
of  8  and  ♦,  which  are  taken  at  5°  intervals.  In  terms  of  the  standard 


J 


stereographic  projection,  0  corresponds  to  the  longitude  and  p  to  the 
latitude;  i.e. ,  ip  =  0  indicates  the  equator  and  \p  -  ±  90°,  the  North 
and  South  poles  of  a  Whiff  net.  Fig.  13  presents  the  mapping  of  E  for 
the  case  when  the  e*j  of  Eq.  (10)  is  used.  The  values  of  k,  E^y/yVy 
and  E^/yV,  are  also  given,  and  k  is  zero  for  Fig.  13a  and  is  0.1 
for  Fig.  13b,  respectively.  Similar  mappings  for  the  e^j  of  Eq.  (11) 
with  k  =  0,  0.1,  0.4  and  0.8,  respectively,  are  shown  in  Figs.  14a 
through  14d.  More  detailed  mappings  are  shown  in  Figs.  15a  and  15b, 
which  present  the  normalized  values  of  E  for  the  range  of 
-10°  £  0  S  20°  and  1^1  ±  10°  at  1°  intervals.  Hie  magnitude  of  E  is 
given  by  an  integer  between  0  and  1000,  where  the  latter  corresponds  to 
the  E  value  at  0  =  -10°  and  ip  =  ±  10°. 

From  Figs.  13  and  14,  it  can  be  seen  that  the  strain  energy  minima 
appear  only  at  the  locations  predicted  in  Figs.  10  and  11  with  p  ~  0. 
However,  the  low  energy  orientation  spreads  over  a  range  of  0  and  ^ 
from  the  minimum  energy  orientation.  This  is  best  illustrated  in 
Figs.  15a  and  b.  E,^  is  located  at  0  =  9°  and  ip  =  0,  but  within  the 
region  of  A0  =  ±  2°  and  =  ±  4° ,  the  value  of  E  increases  over  I^.n 
less  than  1%  of  E  of  a  sphere  with  the  identical  e*j  and  Vj.  Hie 
corresponding  ranges  of  A0  and  tp  for  other  values  of  k  and  e^  are 
listed  in  Table  II.  For  the  two  choices  of  e*j  employed,  the  low 
energy  orientation  far  thin  discs  (k  $  0.1)  is  located  over  a  wider 
range  of  p  than  that  of  0. 


jfiic 


S.  Discussion 


5.1  Effects  of  Rotational  Components 


When  the  transformation  strain  et..  is  given,  Eshelby^10^  showed 
that  the  uniform  rotation  w?.  =  i  (u?  .  -  u?  • )  in  the  inclusion  is 

1J  *  J 


written  as 


wij  =  nijki 


The  only  non-zero  components  of  see  nl212  1I2323,  ni313  i-11  a 

coordinate  system  whose  axes  coincide  the  principal  axes  of  the 
ellipsoidal  inclusion.  These  are  given  in  terms  of  elliptic  intergnals 


by,  e.g.. 


"1212  =  "n2112  =  (Ib  "  V/8ir  * 


Since  the  rotational  components  of  the  transformation  distortion  must 
be  subetracted  in  order  to  apply  the  Eshelfcy  theory,  the  uniform  rotation 
in  the  inclusion  from  the  original  lattice,  ,  is  given  by 
I  c  t 

o»ij  =  -  <*>£j  .  (16  J 

This  rotational  term  most  be  considered  in  comparing  the  results  of  the 
previous  section  to  the  crystallographic  orientations  determined  by 
experiment. 

In  both  the  e^  employed  in  the  above  calculations  of  E,  the  only 

c  t  t 

non-zero  contribution  to  arises  from  e£2  *  e21  =  0.088,  and  is 

«12  =  -w21  =  2  ej2  ni212  -  ej2  . 


(17) 


18 


The  angle  of  rotation  <j>3  about  axis  becomes 
<t>3  =  tan-1  C— <^2^  * 

The  numerical  results  are  given  in  Table  III,  which  shows  that 
the  rotational  correction  amounts  to  only  1.05°  for  k  =  0.1,  and 
vanishes  for  k  =  0. 


TABLE  III 

Rotational  Corrections  (for  e^2  =  0.088) 


k 

ni212 

CM 

<|>3  (degrees) 

0 

0.5000 

0.0 

0.0 

0.1 

0.3956 

-1.837  x  10~2 

1.05 

0.4 

0.1912 

-5.436  x  10"2 

3.11 

0.8 

0.0458 

-7.994  x  10-2 

4.57 

5.2  Habit  Planes  and  Shape  Deformation 

The  minimum  strain  energy  configuration  of  a  spheroidal  inclusion 
corresponds  to  the  most  likely  geometry  of  the  new  phase,  when  the 
interface  energy  contribution  to  the  nucleation  process  can  be  ignored. 
If  the  interface  energy  contributes  significantly,  its  effect  must  be 
evaluated  as  was  done  in  Section  4.1,  However,  the  orientation  depen¬ 
dence  of  the  interface  energy  is  generally  unknown.  If  the  interface 
energy  is  assumed  to  be  independent  of  orientation ,  the  most  favored 
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orientation  of  the  new  phase  is  essentially  governed  by  the  strain 
energy  except  for  the  fact  that  an  optimum  value  of  the  aspect  ratio 
exists  where  the  sum  of  the  strain  energy  and  interface  energy  contri¬ 
butions  is  minimized. 

Considering  the  -transformation  strain  of  Eq.  (11)  that  corresponds 
to  the  c/a-ratio  of  1.586  for  the  product  phase,  the  strain  energy 
minimum  exists,  when  the  broad  face  of  the  spheroidal  inclusion  (k  =  0) 
is  parallel  to  (9  9  13)g.  When  the  aspect  ratio  is  increased  to  0.1,  the 
minimum  strain  energy  orietation  is  shifted  by  4°  to  (9  9  14 )g  including 
the  effect  of  <j>3-  These  orientations  are  shown  in  Fig.  16a  together  with 

the  ranges  of  orientations  where  E  increases  over  E_.  less  than  1%  of 

mm 

T 

E  of  a  sphere  with  identical  e^  and  V^.  These  ranges  are  shown  by 
a  hatched  area  (k  =  0)  and  by  a  shaded  area  (k  =  0),  respectively.  The 
present  result  of  (9  9  13  )g  habit  plane  is  almost  identical  to  experi¬ 
mentally  determined  habit  planes  obtained  by  Newkirk  and  Geisler^15^ 

(16  ) 

and  by  Liu.  The  observed  iabit  plane  orientations  are  indicated 

in  Fig.  16b  by  a  filled  triangle  [the  (8  8  11  )g  orientation]  and  by  a 

hatched  area,  respectively.  The  range  of  observed  habit  plane  orientations 

(17) 

shown  by  a  shaded  area  is  due  to  Gaunt  and  Christian,  whereas  a  double 

(18  ) 

circle  represents  the  result  of  Williams  et  al.  [the  (8  9  12)g  orientation] 

(19) 

and  a  filled  square  that  of  van  Ginneken,  [the  (569)g  orientation]. 

All  the  observed  habit  plane  orientations  are  within  a  few  degrees  of 
the  predicted  low  energy  regions,  indicating  an  excellent  agreement  between 
theory  and  experiment.  It  is  also  significant  that  the  predicted  orienta¬ 
tions  of  low  strain  energy  are  not  confined  to  a  narrow  zone.  This  is 
again  in  acoord  with  experiment,  since  it  has  been  recognized  that  the 
observed  scatter  of  habit  planes  is  beyond  the  limit  of  experimental 


error. 
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The  predicted  habit  plane  of  the  Bowles-Mackenzie  analysis  for 
pure  titanium  is  (8  9  12)g^20^  and  agrees  with,  the  observed  one  reported 
by  Williams  et  al.  However,  the  agreement  hinges  upon  a  suitable 

choice  of  the  adjustable  dilatation  parameter,  which  amounts  to  0.7%.  The 
lattice  invariant  shear  plane  in  the  above  analysis  is  {10ll}^.  Twins 
on  {l0ll}„  in  Ti  martensite  have  been  indeed  observed  by  transmission 
electron  microscopy,  but  their  thickness  is  of  such  magnitude  that 

it  is  difficult  to  interpret  them  as  the  lattice  invariant  shear. 

The  Bowles-Mackenzie  prediction  of  the  habit  plane  is  rather  close  to 
the  predicted  low  energy  orientations  of  the  present  analysis.  This  is 
perhaps  not  surprising  since  it  is  expected  that  the  invariant  plane 
strain  condition  postulated  by  the  Bowles-Mackenzie  analysis  should 
approximate  the  strain  energy  minimum  condition. 

The  total  shape  distortion  in  the  new  phase,  8^,  is  given  by  the  sum 
of  e?.  and  d •  in  the  present  analysis.  For  the  above  case  where  e. •  of 

1J  1J 

T 

Eq.  (11)  is  used,  becomes 


(20a) 


for  the  case  of  k  =  0  and 


(20b) 


for  the  case  of  k  =  0.1. 


^ 'y;.-  ^  "  •  -  titi^j i^^iMq^;!  . "  i 
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The  sum  of  By  and  an  identity  matrix  corresponds  to  the  total  shape 

T 

deformation  F  in  the  crystallographic  theory.  As  a  comparison ,  By 
obtained  from  the  analysis  on  pure  Ti  by  Bowles  and  Mackenzie  is 
given  in  the  same  coordinate  system  by 


b.T 

13 


/ -0.1075 
(  0.0550 
V  -0.0658 


-0.0618 

0.0268 

0.0405 


(21) 


T  . 

It  is  clear  that  By  in  the  Bowles-Mackenzie  analysis  is  not  optimized 

T 

for  the  strain  energy  minimum,  as  the  strain  components  of  By  are  al¬ 
most  universally  greater  than  those  of  e?.  in  the  present  analysis. 

A  comparison  of  the  two  8-4 's  indicates  that,  although  the  net  dilatation 
CBjJ  +  @22  +  differs  only  slightly*  a  large  contraction  along  the 
[Hl]B  and  a  substantial  expansion  along  the  [01i:B  are  present  in  the 
ft.,  for  the  Bowles-Mackenzie  analysis  [Eq.  (21)].  The  corresponding 
dilatation  parameters  in  the  present  analysis  are  zero.  The  origin  of 
such  a  difference  can  be  traced  to  the  processes  of  lattice  transformation 
in  the  crystallographic  theory,  where  the  invariant  plane  strain  condition 
is  accomplished  by  adjusting  the  dilatational  strain  components.  This 
approach  is  in  a  sharp  contrast  to  the  present  one  where  the  transformation 
strain  is  chosen  by  minimizing  atomic  movements  involved.  It  is  also 
evident  that  the  dislocation  concept  of  shear  transformation  is  completely 
absent  in  the  lattice  transformation  stage  of  the  crystallographic  theory. 
It  is  invoked  only  during  the  lattice  invariant  shearing.  This  is  why  it 
is  so  difficult  to  reconcile  the  crystallographic  theory  with  the  atomistic 
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or  dislocation  concept  of  martensitic  transformation.  The  former 
achieves  the  lattice  transformation  via  pure  strain  (e.g.,  the  Bain 
distortion),  whereas  the  latter  contends  the  major  changes  to  be  due  to 
the  motion  of  transformation  dislocations.  In  the  present  theory,  e^ 
arises  from  the  atomistic  concept,  which  is  therefore  an  integgral  part 
of  the  theory. 

5.3  Needle-  and  Lath-Type  Precipitates  of  a-Ti 

When  a  concentrated  solution  of  Ti  with  Mo,  V  and  other  so-called 
3-stabilizers  is  solution -treated  and  aged  in  a  a  +  3  phase  region,  needle- 
and  lath-type  precipitates  of  a-Ti  are  found  in  the  a-Ti  matrix.  The  long 
axis  of  these  precipitates  is  aligned  along  <110>E,  and  the  Burgers  rela¬ 
tionship  is  satisfied  in  most  of  the  instances  investigated/11^  When 
the  aging  was  performed  at  relatively  low  temperatures  (430-480°C  for 
a  Ti-18Mo  alloy  and  480-540°C  for  Beta-Ill  Ti  alloy),  Rosales^1^  also 
noted  that  the  alloy  partitioning  considerably  lags  behind  the  3  to  a 
transformation  and  suggested  a  shear  transformation  mechanism.  Since  the 
lattice  correspondence  is  identical,  the  present  analysis  also  applies  to 
the  a-precipitation .  The  results  presented  in  Section  4.1  indicate  that, 
while  a  disc-shaped  inclusion  has  the  lowest  strain  energy,  a  needle- 
shaped  inclusion  along  <110>g  are  the  preferred  geometry.  It  is  not  known 
what  causes  the  needle-  and  lath-type  precipitates  to  be  most  stable,  but 
it  is  naturally  expected  that  the  diffusion  of  solute  atoms,  the  nature  of 
nucleation  sites  such  as  grain  boundaries,  the  structure  of  a- 3  phase 
boundaries  and  precipitate  growth  mechanism  have  as  much  influence  as  the 
strain  energy.  The  k-dependence  of  E  as  shown  in  Figs.  6-9,  however, 


presents  one  possible  explanation.  If  the  nuclei  of  o-precipitates 

(21) 

have  semispherical  geometry  as  proposed  by  Rosales,  the  growth 
of  <110  needles  is  preferred  over  that  of  <112>g-discs  according  to 
the  present  results  (cf.  Fig.  8).  This  is  due  to  a  higher  value  of 
-dE/dk  for  the  needle  geometry  in  comparison  to  that  for  the  disc 
geometry  when  k  is  close  to  unity;  that  is,  the  rate  of  reduction  of  E 
is  greater  for  the  growth  of  < 110>g-needles . 

5.4  General  Discussion 

The  theory  described  in  this  paper  represents  a  significant  de¬ 
parture  from  the  crystallographic  theory  or  from  the  interface  dislocation 
(22  23) 

models.  ’  It  combines  the  atomistic  transformation  mechanism  based 
on  dislocation  movement  and  the  crystallographic  features  of  transforma¬ 
tion  through  the  suitable  selection  of  the  transformation  strain  tensor. 

The  theory  can  be  extended  to  othernartensitic  transformations ,  some  of 
which  are  being  investigated  by  the  present  authors. 

The  present  theory  can  be  refined  by  considering  the  elastic  constants 
of  an  inclusion,  which  differ  from  those  of  the  matrix.  Ihis  refinement 
can  be  done  within  the  framework  of  the  Eshelby  theory.  The  next 
extension  of  the  theory  will  involve  the  consideration  of  elastic  aniso¬ 
tropy.  The  general  theory  for  anisotropic  inclusion  problems  has  been 
developed  by  Kinoshita  and  Mura/ 24 ^  Although  a  few  specific  solutions 
for  a  spheroidal  inclusion  have  been  reduced  to  line  integrals/25^  these 
are  not  inmediately  applicable  to  martensitic  trans f orrot ion  problems. 

This  is  due  to  the  symmetry  requirement  for  the  elastic  constants.  Only 
the  spheroidal  inclusions  whose  major  axes  coincide  writh  the  cube  axes  can 
be  treated  by  the  solutions  available. 
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An  entirely  different  approach  to  the  transformation  problem  has 
been  developed  by  Khachaturyan.  '  The  theory  utilizes  the  computation 
of  elastic  strain  energy  in  the  reciprocal  space.  A  general  solution  is 
obtained  for  an  inclusion  problem  in  the  anisotropic  elasticity  and  a 
specific  solution  is  given  for  a  tetragonal  dilatation  in  the  cubic 
symmetry.  A  similar  approach  has  been  successful  in  the  analysis  of 

f  OC  0*7  } 

the  omega  transformation  in  Ti  and  Zr  alloys.  *  '  However,  the 

application  of  this  theory  to  martensitic  transformation  problems 
requires  a  substantial  effort  in  the  futuare. 

6 .  Conclusions 

1.  The  total  strain  energy  in  matrix  and  a  spheroidal  inclusion 
is  determined  by  using  the  Eshelby  theory.  Approximate  formulae  and  a 
parametric  expression  for  E  are  obtained  by  evaluating  elliptic  integrals. 
Numerical  results  for  an  oblate  spheroid  are  included  for  the  most  general 
form  of  the  transformation  strain,  e^j. 

2.  The  atomistic  transformation  mechanism  is  combined  with  the 
Eshelby  theory  to  provide  a  new  approach  to  the  analysis  of  martensitic 
transformation.  The  correspondence  and  distinctions  between  this  new 
theory  and  the  phenomenological  crystallographic  theory  are  clarified. 

3.  Specific  e^j  for  the  bcc-to-hcp  transformation  in  Ti  and  its 
alloys  is  obtained.  The  habit  plane  is  predicted  on  the  basis  of  the 
strain  energy  (and  interface  energy)  minimization  principle.  The  results 
are  in  excellent  agreement  with  experiment. 
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FIGURE  CAPTIONS 


Schematic  representation  of  the  Eshelby  theory and  the 
corresponding  processes  in  the  crystallographic  theory. 

Coordinate  axes  employed  and  an  oblate  spheroid  with  the  broad 
face  //(5!l)g. 

Strain  energy  coefficients  for  an  oblate  spheroidal  inclusion. 


The  rotation  of  the  coordinate  axes  about  the  Xg  axis. 

Atomic  arrangements  on  {110}g,  {110}g  after  shearing  and  (0001)^  . 
Burgers  relationship  is  maintained  in  the  atomic  movement. 

Strain  energy  E^  due  to  a  shear  component  e^  against  k. 

Strain  energy  E ^  due  to  a  dilatation  component  e^  against  k. 

"t 

Total  strain  energy  E  against  k  e. .  of  Eq.  (10)  was  used  and  E  is 

t  ” 

expressed  in  terms  of  . 


The  sum  of  E  and  the  interfacial  energy,  Ep, against  k. 


The  variation  of  E  due  to  rotation  about  the  Xg  axis,  using 
of  Eq.  (10).  k  values  are  given  in  the  figure,  and  0=0 
corresponds  to  [5Il]-disc  and  0  =  90°  to  [lll]-disc. 


The  variation  of  E  due  to  rotation  about  the 
of  Eq.  (11). 


t 

x~  axis,  using  e.. 
3  lj 


The  k-dependenoe  of  E  .  .  et .  used  is  indicated. 

min  ij 

The  strain  energy  napping  with  respect  to  0  and  e^j  is  given 
by  Eq.  (10).  locations  are  given  by  (*)  and  Erni-n  by  (+). 

(a)  k  =  0  (b)  k  =  0.1. 


The  strain  energy  mapping  using  e^  of  Eq.  (11).  (a)  k  =  0, 

(b)  k  =  0.1  (c)  k  =  0.4  (d)  k  =  0.8. 


Details  of  Figs.  14a  and  b,  in  the  vicinity  of  0  =  =  0. 

(a)  Habit  planes  predicted  by  the  present  theory 

(b)  Experimentally  determined  habit  planes. 


energy  E,  due  to  a  shear  component 
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lue  to  a  dilatation  component  against  k. 
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Total  strain  energy  E  against  k 
expressed  in  terms  of  e*„. 


variation  of  E  due  to  rotation  about  t 
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Fig.  13  The  strain  energy  napping  with  respect  to  9  and  \|>.  e^  is  given 
by  Eq.  (10).  locations  are  given  by  (*)  and  by  (+). 

(a)  k  =  0  (b)  k  =  0.1. 
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